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Abstract 

Simulations in compact U(l) lattice gauge theory in 4D show now beyond any reason- 
able doubts that the phase transition separating the Coulomb from the confined phase 
is of first order, albeit a very weak one. This settles the issue from the numerical side. 
On the analytical side, it was suggested some time ago, based on the qualitative anal- 
ogy between the phase diagram of such a model and the one of scalar QED obtained by 
soft breaking the N = 2 Seiberg-Witten model down to N — 0, that the phase tran- 
sition should be of second order. In this work we take a fresh look at this issue and 
show that a proper implementation of the Seiberg-Witten model below the supersym- 
metry breaking scale requires considering some new radiative corrections. Through the 
Coleman- Weinberg mechanism this turns the second order transition into a weakly first 
order one, in agreement with the numerical results. We comment on several other aspects 
of this continuum model. 
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1 Introduction 



£7(1) pure gauge theory in four dimensions is an interacting theory when formulated on 
the lattice. In fact there are two physically distinct phases: one characterized by massless 
photons and long-range interactions (the regime we would associate to continuum QED) and 
a strongly interacting region, where electrical charges are confined and monopoles (that have 
a finite energy thanks to the finite value of the lattice spacing) proliferate. It is a priori 
unclear which "continuum" theory is adequate to describe this regime. 

A long-standing controversy regarding the order of the phase transition separating both 
regimes existed for some time. Traditionally the phase transition was believed to be of first 
order [Q, but this belief was questioned in a series of papers where some evidence was 
provided to suggest a second order transition. These authors believed that the apparent first 
order nature of the phase transition was due to some specific configurations that on a toroidal 
topology would lead to long-lived metastable states and these could be mistaken as footprints 
of a first order transition. To avoid metastability the authors performed a simulation using 
spherical topology. The same authors also performed a number of interesting measurements 
regarding the spectrum of the theory [3]. 

Soon afterwards, it was realized that, when going to larger lattices, a signal for a first 
order transition was seen also on spherical topologies @J, thus concluding that the transition 
was first order even if a very weak one. This seems to be the commonly accepted lore [5]. 

Clearly in all the above discussion an analytical approach is lacking. It would be highly 
desirable to gain some understanding of the mechanism of monopole condensation in this 
specific problem, but until very recently it was unclear how to treat these topological ex- 
citations. Furthermore, the fact that the transition was believed to be of first order would 
suggest that indeed there was no continuum theory associated to it. The suggestion that the 
transition could be of second order, combined with the results on N = 2 and N = 1 gauge 
theories obtained by Seiberg and Witten [H] led a group of authors to propose a continuum 
model for it and, in fact, for the strongly coupled regime of compact lattice £7(1). 

We shortly summarize the scenario and results presented in [7]. We shall refer in what 
follows exclusively to the continuum model. The starting point is N = 2 SU(2) Yang-Mills 
theory, whose action is 



'bare 
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Its scalar potential has a flat direction and in order to define the theory a complex parameter, 



the vacuum expectation value of the scalar field, has to be fixed, ((f)) = har 3 . 



At scales /x 2 <C u = ha 2 the charged fields decouple and the model is effectively described 
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by a N = 2, U(l) gauge theory 
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where F is a known function of A (the chiral superfield in the N = 2 multiplet of the photon) 
and the dynamically generated scale A, W a is the Abelian field strength (for details and 
notation we refer to [SJ |H| ) . For scales fj, 2 > u the charged fields are not decoupled and they 
contribute to the running. At /i 2 ~ u, when the charged fields decouple, the effective coupling 
constant, which is related to the function F by 47r/<7 2 = Im^^-, freezes at the value g = g{u) 
since below /j, 2 = u no fields can contribute to its running. If u 3> A, this value can be safely 
computed within perturbation theory by using the beta function of N = 2, SU(2) Yang-Mills 
theory j^]. The value g = g(u) is thus the appropriate one to use in the effective potential at 
p 2 0. 

As u decreases and approaches A, perturbation theory becomes unreliable. The appro- 
priate variables to describe the theory are not g and a, but the dual JU] variables go and 
ao, defined through the relations 

OF Att T d 2 F , 
aD = Ta -J D = Im dal = hl (3) 

Holomorphy and monodromy dictate that gr> runs around go ~ in a way that shows that a 
massless hypermultiplet is present in the spectrum: it corresponds to a magnetic monopole. 
The appropriate value of the coupling constant for long-distance physics (p 2 — > 0) in this 
regime is goi^D)- This is exactly for the same reasons that were described above phrased 
this time in term of dual variables. 

The appropriate effective theory description for small values of od is thus provided by [S] 
S = Sp =2 + Sm, where 

S#= 2 = S SW (A^A D ,V^V D ), 
S M = I d 4 x J d A 9 ^M*e 2VD M + M*e- 2VD M^J + J d 2 6V2A D MM + h.c. (4) 

and where M and M describe the N = 2 monopole hypermultiplet. Corrections to the above 
effective action will be of order p 2 / A 2 . 

So far this seems to have little to do with the compact lattice U{\) theory, except for the 
fact that the theory has a manifest U(l) symmetry and that massless monopoles appear in 
a particular point of the moduli space of the theory. Clearly we need two more ingredients: 
we need to enlarge the region where light monopoles appear in order to get a finite density 
of monopoles (m) 7^ and, eventually, to break N = 2 down to N = 0. 

In [7] supersymmetry is broken in two steps. First one adds a coupling between the chiral 
multiplet $ (a member of the N = 2 hypermultiplet) and a N = 1 chiral superfield z; as a 
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consequence the model has now residual N = 1 only. The part of the action containing the 
new super field is: 



(5) 



It introduces two free parameters I and w. This term has the net effect of enlarging the 
monopole condensation region. 

The limit / —* oo used in [7j to study the vacuum structure allows to integrate the z super 
field out using its classical equations of motion. Unlike the original mechanism of Seiberg 
and Witten [H] the breaking introduced by the term S z is a hard one 1 and introduces some 
non-trivial quantum corrections. The obvious main effect of adding S z to the action is to lift 
the degeneracy of the vacuum so u = ^a 2 is no longer a free parameter; in fact when taking 
the limit I — > oo it will be kept close to w. Again, if u (now w) is large, perturbation theory 
is valid and a semiclassical calculation makes sense. The coupling constant now will run with 
the perturbative beta functions of N = 1 (extended with the new matter multiplet z in the 
case of I — ► 0). However at some point, exactly as for N = 2, the perturbative procedure will 
break down. 

Then to go from N = 1 to N = and make contact with the real world, the technique 
discussed in [H] was used. This consists in coupling to the original N = 2 superfield a further 
N = 2 superfield (spurion). The action becomes 



~,N=0 1 t 

>sw = —Im 



4n 



1**1 



1 J 



i,j = 0,l. 

(6) 

The physical fields are understood to be those labeled '1'. Breaking to N = can be achieved 
by giving non-zero values either to the D or F terms of the spurion (labeled '0'). In [7], the 
choice of the auxiliary fields to break to N = is Dq ^ 0, -Fo = for the reasons mentioned 
there. The starting point is thus finally S z + Sg^° 

We are now in a position to make contact with compact lattice U(l). Let us first sum- 
marize the parameters that we have at our disposal. If we forget about the parameter I 
(we simply assume that it is large enough to enforce the constraint u ~ w), we have two 
dimensional parameters, namely w and Dq and a dimensionless one, g, that can be traded 
by A. If we are able to place ourselves in a region where the characteristic momenta is much 
below both w and Dq it is clear that we will be dealing with an effective N = 0, U(l) gauge 
theory. We still have some freedom in adjusting the relative values of w and Dq; if by doing 
so we are able to trigger monopole condensation we shall be describing the continuum version 
of the confinement-Coulomb phase transition seen in compact lattice U(l). If, in addition, 
1 It does not appear possible to enlarge the monopole condensation region with a soft breaking. 
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we manage to do all the previous manipulations in a controlled manner some quantitative 
predictions can be made. It is our purpose to convince the reader that this is the case. 



2 The necessity of new quantum corrections 

When constructing supersymmetric effective actions, the non-renormalization theoremsjllj 
determine which perturbative and non perturbative quantum corrections are possible. The 
class of possible counter terms is greatly restricted when supersymmetry is broken by soft 
terms, like in the original proposal to go from N = 2 to N = 1 of Seiberg and Wittenj6j . 

In [S] iV = 2 is broken down to N = by the spurion mechanism allowing for a deter- 
mination of the effective potential of the softly broken theory. (In fact, the first term of an 
expansion in powers of -A around the unbroken solution is found. ) 

The scenario studied in [7j is somewhat different. There the N = 2 supersymmetry of the 
Seiberg- Witten model is broken down to N = by a combination of the spurion mechanism 
we have just discussed and a N = 1 hard breaking term. This one — unlike the breaking of 
Seiberg- Witten or the spurion mechanism — changes the monodromy of the original N = 2 
theory non-trivially. 

This issue was discussed in great detail in 0. By considering the U(V)a x U(1)r charges, 
it is possible to see (modulo some highly plausible assumptions) that the requirement of 
N = 1 suffices to constrain the effective action in a way that is sufficient for our purposes. 

Since we are concerned about the monopole condensation mechanism in compact U(l) 
we must be in a region of moduli space close to A. This requires, as explained in jS], the use 
of the dual version of the Seiberg- Witten effective action Ssw completed with a term for the 
monopole and the dualization of S z , 

S = S%=° + S M + S?, (7) 

S%=° = S^°(A - A D , V - V D ), (8) 

5f = J d 4 x(^J d 4 9K(z,J) J d 2 9lzS.w-U(A D ) + ^f(^,^\+h.c?j (9) 

This involves an unknown function as well as an undetermined Kahler potential (originated 
by the quantum effects of the N = 1 term arising when constructing the effective action [7]) 
which are not calculable by symmetry considerations . We shall assume that the net effect 
of this term when computing the scalar potential is to adjust the value of a£>; namely we 
exchange w by old as a free adjustable parameter. 
2 As usual, terms of order p 2 /A 2 are neglected. 



where 
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This, in our understanding, does not exhaust the quantum corrections we have to include. 
At this point we depart from the analysis done in |7J. Let us see why the analysis presented 
in [7] is incomplete. 

We have been quite careful about the actual meaning of the coupling constants appearing 
in the effective action. They are constants renormalized at the scale ap, when the dual 
variables are used. This is fine as long as the coupling constants and parameters in the 
effective Lagrangian do not run from the scale ar> to p 2 = where the effective potential is 
defined. This is the case in the N = 2 theory, or even if the theory is broken to N = 1 or 
N = when the breaking occurs at a scale well below ap. In the present case supersymmetry 
is severely broken (i.e. p 2 -C Dq ~ ar>) and the running needs to be taken into account. This 
makes the effective potential subject to radiative corrections. 

We argue, in fact, that starting from energies around Dq new quantum effects (the stan- 
dard ones for non supersymmetric QFT) have to be included. To be more precise, the 
hierarchy of scales we have is 

< p 2 < b 01 D ~ a D < A. (10) 

Starting from the left, the first inequality is dictated by the interest in studying the non 
supersymmetric regime, whereas the second allows us to use the results of L 8 . In fact we will 
use the analytical solution of Seiberg-Witten and its generalization to iV = in (J7J) as the 
Landau-Ginzburg (tree level) approximation to the complete effective action . From there 
we still need to scale down to zero momenta to read out the effective potential. By doing 
so we shall show how the agreement between the conjecture in [7j and the lattice results 
is recovered. In fact, the Coleman Weinberg mechanism takes place: radiative corrections 
transform a second order phase transition into a first order one ^2] • 

3 The classical vacuum structure and mass spectrum 

After elimination of the auxiliary fields and the z superfield (through the limit I — > oo) we 
are left with a scalar potential that, up to constant terms, reads 

V = — — (rh)m — m)m) + - — |a£)| 2 {ram} + mfhf\ + -^-Dq (mm^ — rh^rh) (11) 
2&ii V / b\i \ / Oxi v / 

In this expression 601 and »oo have to be understood as functions of the point of the moduli 

(their analytical expressions can be found in |S]) and b\\ is the gauge coupling (derived 

in |S]). The different supersymmetry breaking terms do not alter the dependence of bij on 

on as shown in [7] above the scale -Do- As discussed in the previous section, we take the 

effective potential (|ll(l as defined at the scale old- In ((TTj) the quantities a£> and Do are free 

parameters. 
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After rotating the fields using the U(l) ® U(l) rigid symmetry of the above potential, it 
depends only on two real fields m and rh. The potential can have, depending on the values 

P9i i 1 2 i 1 2 

of the parameters, three different minima 7 : m = rh = when — 2 |od| < < 2 |od| , 

or m = and rh 2 = — (2b\\ (z 2 + z|) — 601-D0) when &01A) > 2 |a£)| 2 , or rh = and m 2 = 

1 1 2 1 1 1 2 

2 dj) + 6oiZ?o when < — 2 |a£)| . 



Here we will focus on the last case (The second and the third one are interchangeable, 
depending on the sign of Dq). In order to study the critical region we introduce as control 
parameter the combination 

a = -(2\a D \ 2 + b 01 D ) . (12) 



M 2 = 

01 



In fact, the mass of the monopole field (responsible for the critical behavior) is given by 

J L(2\a D \ 2 + b 01 D ) (13) 

and is thus proportional to a. We demand (to stay close to the transition) a — > 0, while 
both Do and cld are understood to be large compared to the physical scales (to decouple 
the supersymmetric modes), but small compared to A (in order to trust the Seiberg-Witten 
effective action). 

In the region we are considering the only light scalar field is the monopole. The rh field 
has a large mass 

M?=^< (14) 

Oil 

and it can be safely eliminated from the effective potential which is just 

V = HiT ( m ? + m s) 2 + T~ ( m 2 + m i) ( 15 ) 
lb\\ On 

where we have restored the imaginary part of the monopole field mi- This is of course an 
ordinary A(/> 4 potential. 

As for the fermions, we have the following non diagonal mass matrix 

/ an ia \ 





2b 
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ClD 










a 





a 





— ia 





1 

16vr 



J_ ( 

167r I d 2 a\daQ 



Mr. \ t/~a I (/tl n 1 ^0 



(16) 



where the rows and columns are ordered as ip m , ipfh, ipi ^1 (see [S] for notation). This matrix 
can be diagonalized to obtain two physical fermions with masses close to and two with 
masses closed to y^U ( gP^ ~da ) groups are extremely heavy and can be dropped 

from the effective action using their equation of motion. 
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4 One-Loop effective potential and phase diagram 



Using the results of Seiberg and Witten as 'boundary conditions' at the scale ao ~ &01A) 
we are now in a position to run the effective potential down from fx = ao to /J, = 0. The 
only light fields that are left are as expected the dual photon and the monopole field. This is 
however not a generic potential: some peculiar relations between self and electric couplings 
are inherited from supersymmetry. But of course these relations shall not be preserved as 
supersymmetry is broken below Dq ~ ao- 

The full Lagrangian, counter terms included, will be 

C = (d»m 1 + -±=m 2 AA + ( d„m 2 - A„m^j + ± ((d 2 g^ - A„A V ) 



m l +m 2 ) — - — [m 2 + m l ) (17) 



2b n " 1 " b n 
+ 2bJ~ ^ + m ^ 2 + 26be (^™i^™2) - 25b e (S^mi) 
+—=4 (m 2 . + ml) A^A M - — — (m| + m{) . 



/Oil °11 

We shall work in the Landau gauge. 



The one loop correction to the effective potential in MS produces the following effective 
potential 

1 A a 2 l/3 2 a \ 2 / 3m 2 + a 3 N 

V = — — m + - — m H ^ - — m + - — log ^ 

26n &ii (4vr) 2 \hi bli/ V ^ 2 . 



1 2 a \ 2 / m 2 + a 3\ 3 m 4 / m? 5 
-m + — log ^ + ——--pr log 



(4^) 2 V^ii buy V V 2 2 J (4i) 2 &?i\ A* 2 6, 
in agreement with for instance. We have used again the residual U(l) rigid symmetry 
to get rid of m 2 and we have replaced m\ — > m. Here fj? is the subtraction scale. 
From the one loop calculations we also extract the beta functions 

3 

1 13 1 / 2 \ 2 

A " = (S?^ *- = 55?(5r) ■ <19> 

which are of course different. 

With these results what follows is quite predictable. We are interested in the phase 
transition from the Coulomb to the confined phase. The order parameter is m. A value of 
m different from zero reveals monopole condensation and the presence of the dual Meissner 
effect. 

By changing a we move from a phase with m ^ to one with m = . If we tune a to 
have zero renormalized mass (to stay on the critical surface) the presence of a minimum in 
the potential outside from the origin implies a first order phase transition [12 . 
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On the critical surface, there are in the effective theory two distinct coupling constants 
(the gauge coupling 6 e and the quartic self coupling for the scalar field £>xi) that have the same 
bare value , as dictated by supersymmetry. We call this bare value bf± . In spite of this, they 
are different below cld as they run with different beta functions (a naive inspection of the 
bare potential would lead to the erroneous conclusion that there is just one coupling). The 
existence of two beta functions, as pointed out by Yamagishi is the essential ingredient 
that allows, in scalar QED, (as in other models ^S]) radiative corrections to induce a first 
order transition ^2]- Following the analysis in ^1] we obtain the one loop RG improved 
effective potential 



v b n (t) 4 

V = ; — m exp 

4! 



dt'- 



7(0 



(20) 



l- 7 (tf). 

where t = ln(m/[j,) and j(t) = j(bu(t),b e (t)) with b e (t),bn(t) solution of the differential 
equations 

db u (t) _ p bll (b n ,b e ) db e {t) _ f3 be (b n ,b e ) 

dt 1- 7(611, b e )' dt 1- 7(611, b e ) 

with initial conditions 

6u(0) = 6ff 6 e (0)=6fr (22) 

The existence of a minimum away from the origin implies that the trajectories (|21jl in the 
(6n, 6 e ) plane should cross the line 

4&n + f3 bll = 0. (23) 

The plot of the renormalization-group trajectories of the coupling constants (tangent to 
the arrow field in figure ^) shows that the trajectories leaving the bare curve j~ = ^ (upper 
solid line in figure ^) always cross the line ()23j) (lower solid line in figure This implies that 
the transition is of first order in all the plane but the origin. A zoom of the previous plot in 
the perturbative region with explicitly sketched flow is drawn in figure 

The fact that the transition is first order can also be checked by direct inspection of the 
potential (|18|) : that is by plotting (|18|) for typical values of the parameters and looking to the 
position of the minima. 

It is evident that the system passes from the Coulomb phase to the confined one through 
a first order transition (fig. We thus conclude that the transition is indeed of first order. 
However, as befits a radiatively induced first order transition, this is expected to be quite 
weak, exactly as observed in lattice simulations. 



5 Conclusions 

The agreement between the predictions of Seiberg-Witten theory and what is observed in 
lattice simulations of compact U(l) goes beyond the fact that both indicate the existence of 
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a weakly first order transition. As it was emphasized in j7j, the parity and charge-conjugation 
quantum numbers of the lightest one-particle states in the confinement phase (0 ++ and l" 1 ) 
seem to agree with the clearest signal of "gauge-ball" state from lattice simulation 3 . The 
fact that the transition is weekly first order makes the comparison between continuum and 
lattice results meaningful. 

It might seem that, at the end of the day after the successive supersymmetry breakings, 
very little is left of the bare Lagrangian with N = 2, SU(2) Yang-Mills gauge invariance 
and that one is left with the Lagrangian of scalar electrodynamics. This perception is indeed 
correct, but it was by no means obvious that this was the appropriate theoretical framework 
to describe the confining-Coulomb transition in compact lattice U(l). It was not even obvious 
that the appropriate theory could be written in terms of local variables. The good theoretical 
control provided by the Seiberg-Witten model allows us to be on firm ground. 

The supersymmetry breaking is well under control. The problem we are interested is 
considerably easier than trying to get predictions for, say, SU(2) gauge theory, because we 
are in a kinematical regime where it is consistent to have Dq <C A. As far as we can see, the 
theoretical prediction is quite robust. 

The results obtained here modify the conclusions of j7j. Some radiative corrections, 
necessary below the supersymmetry breaking scale, were overlooked in that work. Their 
consideration reconciles the qualitative agreement between the phase diagram of the U(l) 
theory obtained after breaking the supersymmetry of the Seiberg-Witten theory and the phase 
diagram of compact U(l) lattice theory in four dimensions. This also supports the idea that 
confinement in four dimensions is due to monopole condensation as already demonstrated 
in the three-dimensional case |16j . We are currently considering the application of the same 
techniques to other models. 
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Figure 1: Flux of the couplings along the renormalization group trajectories, b\\ versus 
b e . Starting with bare conditions on the upper line the flux always crosses the lower line 
representing 1231 This is the footprint of the first order transition 
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Figure 2: Zoom of the previous plot in the perturbative region (small b e ). The continuous 
line sketches the flux of the couplings along the renormalization group trajectories in the 
bn,b e plane for some typical initial conditions 
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Figure 3: Dynamical evolution: from Coulomb to confined through a first order transition. 
Potential as function of m and a. 
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Figure 4: First order phase transition: the potential as function of m and rh with a = 0.008 
shows the coexistence of the Coulomb and confined phases. 
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